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Abstract 



As it follows from the classical analysis, the typical final state of 
the dark energy universe where dominant energy condition is violated 
is finite time, sudden future singularity (Big Rip). For a number of 
r~| ! dark energy universes (including scalar phantom and effective phan- 

tom theories as well as specific quintessence model) we demonstrate 
that quantum effects play the dominant role near Big Rip, driving 
the universe out of future singularity (or, at least, moderate it). As 
\ a consequence, the entropy bounds with quantum corrections become 

well-defined near Big Rip. Similarly, black holes mass loss due to phan- 
tom accretion is not so dramatic as it was expected: masses do not 
vanish to zero due to transient character of phantom evolution stage. 
Some examples of cosmological evolution for negative, time-dependent 
equation of state are also considered with the same conclusions. The 
application of negative entropy (or negative temparature) occurence 
in the phantom thermodynamics is briefly discussed. 
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1 Introduction 



Recent astrophysical data, ranging from WMAP observations to high red- 
shifts surveys of supernovae, indicate that about 70 percent of the total 
energy of the universe is to be attributed to a strange cosmic fluid with neg- 
ative pressure, dark energy. ^From another side, the universe is accelerating 
currently. It is also observed that the equation of state parameter w is close 
to —1, most probably being below —1. (The possibility of time-dependent 
negative w is not excluded too).^ 

The case with w less than —1 is often dubbed as phantom dark energy. 
At the moment, there is no satisfactory theoretical description of phantom 
dark energy (for a number of attempts in this direction, see [1, 2] and for 
recent review of dark energy, see[3]). The easiest current model of phantom 
is motivated by quintessence [4], it is just scalar field with wrong sign for ki- 
netic energy term. Definitely, such theory being instable shows some weird 
properties caused mainly by the violation of dominant energy condition. In- 
deed, the energy density grows with time in the phantom universe so that in 
a finite time such universe ends up in the singularity dubbed as Big Rip [5] 
(see also earlier discussion of finite time singularity in [6]). The related phe- 
nomenon is that all black holes loss their masses to vanish exactly in Big Rip 
[7]. Phantom thermodynamics looks also strange leading to negative entropy 
of the universe [8] (and divergent entropies near Big Rip) or to appearence 
of negative temparatures [9]. If our universe is indeed phantom one, this all 
may call to revision of basic physical principles governing our reality! 

In the present article aiming to discuss the final state of (phantom) dark 
energy universe we show that situation is much less dramatic than it looks 
from the very beginning (even in the absence of consistent phantom theory). 
Indeed, with the growth of phantom energy density the typical energies and 
curvature invariants grow as well. As a result, much before Big Rip the 
quantum effects start to play the dominant role. In a sense, second Quan- 
tum Gravity era begins. The simple account of the quantum effects, in the 
same way as it was proposed in refs.[10, 11], demonstrates that Big Rip sin- 
gularity is moderated or even does not occur at all. As a result, the entropy 
bounds remain to 1)0 meaningful and black holes masses do not vanish to 

^It is quite possible that what looks like dark energy is the manifestation of some 
unknown feature of the gravitational theory which apparantly should be modified. 
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zero. This observation indicates also that phantom stage (if it is reahstic) is 
just transient period in the universe evolution. 

The paper is organized as follows. In the next section we start from the 
finite time, sudden singularity model proposed by Barrow [12] and consider 
its generalization and give its lagrangian description in terms of scalar-tensor 
theory. In fact, in terms of scalar-tensor theory the sudden singularity is 
different from the model of [12]. It is interesting that such model where 
dominant energy condition is also violated is not necessary phantom with 
wrong sign for kinetic term. Then we show that the account of quantum 
effects (using quantum energy density and pressure obtained by integration of 
the conformal anomaly) moderates the finite time singularity or even prevents 
it. The universe presumably ends up in deSitter phase (future infiationary 
era). A similar analysis has been done in [10] but with coefficient of one of 
the terms in conformal anomaly (OR) being equal to zero. In the present 
paper, due to the importance of the corresponding term at high energies 
(as being advocated by Hawking) the account of arbitrary such coefficient 
is made. Moreover, the backreaction from the matter is included. Section 
three is devoted to the study of final stage for scalar phantom universe and 
effective phantom universe [13] produced by higher derivative coupling of 
scalar kinetic energy with curvature. Again, quantum escape of Big Rip 
occurs or, at least, the singularity is moderating permitting the evolution 
after Big Rip time. In section four the entropy bounds near Big Rip are 
studied. Basically, the typical entropies (including the one for negative time- 
dependent equation of state universe) diverge at singularity. The account 
of quantum effects makes the entropies to be finite and the entropy bounds 
to be well-behaved. Section five is devoted to study of more general, time- 
dependent equation of state which may be also effectively phantom. The 
examples where scale factor is accelerating are presented and the occurence 
of Big Rip is mentioned again. Finally, some summary and outlook are given 
in Discussion. In Appendix A the entropy is written for specific model of 
phantom thermodynamics. It becomes negative for positive temperatures, 
and positive if temperatures are negative. In Appendix B the evolution of 
black holes mass in phantom universe is discussed, also in the case when 
equation of state is time-dependent. The same quantum effects which drive 
the final state out of Big Rip significally improve the evolution of black 
holes mass. It may increase or decrease by phantom energy accretion but 
eventually does not vanish. 
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2 Lagrangian description of classical sudden 
future singularity and quantum effects ac- 
count 

In this section, we construct the scalar-tensor theory with specific potential 
which describes classical sudden future singularity. The account of quantum 
effects near to sudden singularity (where future quantum gravity era starts) 
shows that sudden singularity most probably never occurs. 

In [12], it has been shown that even if the strong energy condition 

p>0, p + 3p>0 (1) 

for some kind of (exotic) matter is satisfied, the future singularity ( Big Rip) 
can occur. Note that dominant energy condition is violated [14] for such a 
scenario. Here p and p are the energy density and the pressure of the matter, 
respectively. We now consider the spatially-flat FRW metric 

ds^ = -dt^ + a{tf {dx'f . (2) 

?=1,2,3 

Following [12], the scale factor a{t) is chosen as 

a{t)^A + Bt'i + C{ts-tf . (3) 

Here A > 0, S > 0, g > 0, and > are constants and C — —At~'^. It is 
assumed that t < tg and 2 > n > 1. There is a singularity at t — > is, where 
^ -Foo. Classical FRW equations 

_6_ 2^ 2 d^g ^ p + 3p 

^' K-'adf 6 ' ^' 

show that 

- - > , p ~ > , (5) 



^ /«2 (^ + Btl) - " ' ^ A + Bt 

near the sing ularity t ^ t,. In (4), if = if . Thus, the energy density p is 
finite but the pressure p diverges. Nevertheless, the strong energy condition 
(1) is satisfied since p and p are positive. 
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Let us present the generalization of Barrow's model [12] as it was pre- 
sented in ref.[10]. In this model, the matter has been given implicitly via the 
FRW equations: 

One may assume H has the following form: 

H{t)^H{t) + A!\ts-t\'' . (7) 

Here H is a smooth, differentiablc (infinite number of times, in principle) 
function and A' and ts are constants. Another assumption is that a constant 
a is not a positive integer. Then H{t) has a singularity at t = ts- In case a 
is negative integer, the singularity is pole. Even if a is positive, in case a is 
not an integer, there appears a singularity, that is, if we analytically continue 
t to the region t > ts from the region t < tg, H might become complex and 
double-valued due to the cut which appears when we analytically continue t 
to be a complex number. It is important for us that singularity presents. If 
we can consider the region t > tg, there is no finite-time future singularity. 
When a > 1, one gets 

pr^-pr^-^H{tsf . (8) 

Hence w — ^ — —1, which may correspond to the positive cosmological 
constant. 

The case < ck < 1 corresponds to Barrow's model and when i ~ i^, we 
find 

pr^—H{ts) , — ^\ts-t\ . 9) 

Here the plus sign in ± corresponds to t < tg case and the minus one to 
t > tg. In the following, the upper (lower) sign always corresponds to t < tg 
{t > ts). The parameter of equation of state w is given by 

^2 A'a\ts-tr' 
w ^ ±- ^ . (10) 

Hence, w is positive in two cases: one is A' > and t < tg, which directly 
corresponds to Barrow's model, and another is A' < t > tg- In other cases, 
w is negative. 
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When — 1 < a < 0, the energy density p and the pressure is given by 



6^",, ,|2a 

P=^^\'ts-t\ , p^± —\ts-t\ . (11) 



The parameter of equation of state is 

w^±^\t,-tr-' , (12) 

which diverges at t = ts- Here w is positive when A' > and t < tg or A' < 
and t > ts- The former case corresponds to sudden future singularity even if 
w is positive. The singularity can be regarded as a Big Rip. (For the recent 
comparison of phantom Big Rip with above type of it, see [15].) The account 
of quantum effects leads to the escape from the future singularity as it was 
shown in refs.[10, 11]. 
The case a = —1 gives 

P^^\ts-t\-\ p^-\ (±2A + 3A'^) % -t\-\ (13) 

which may correspond to the scalar field with exponential potential. The 
parameter w is given by 

= -1 T -| . (14) 

Near t — tg, the universe is expanding if A' > and t < is or A' < and 

t > ts- The former case corresponds to the phantom with w < — 1. In the 
latter case, if 2 > A' > 0, the equation of state describes the usual matter 
with positive w and if A' > 2, the matter may be the quintessence with 
> w > -1. 

If a < — 1, one obtains 



6A 



/2 



P^-P^ -^{ts-tl " , (15) 

which gives w = —1 as for the cosmological constant case. In this case, 
however, there is sharp singularity at t = ts since both of p and p diverge at 
t — ts- This is contrary to the case o; > — 1 as they tend to infinity in the 
limit t ^ ts more rapidly than in case a > —1. 
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Let us construct the Lagrangian (scalar-tensor) model which contains 
sudden future singularity. We start from rather general action of scalar field 
(f) coupled with gravity: 



Then the energy density p and the pressure p are given by 



1 fdchV 1 ' ''^ ^ 



/? (/3 - 1) 00 % - tf-' = /? (/3 - 1) K<l> ' 
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it follows 



(16) 



The scalar equation of motion is 

We are searching for the potential V{(f)), which gives a solution 

H^ho + h (ts -ty , (l)^<f>o (t, - tf (19) 
with constants /iq, /ii, a, and /3. Since 

(P6 



H = hQ + hifJ(t)'l . (20) 
^Prom the FRW equation 



^2 2\dt 

^ (hi + 2h^h^fJr^ + - Y'/'l'^'"^ , (22) 



which gives 

12a 



V i<P) = ^ (hohf.^J-' + hl4>~ ^ <p^-'j -iP-i) P4>U'-^ . (23) 

On the other hand, from the scalar equation of motion: 
one obtains 

--'(^) - -§--f 

Comparing (23) and (25), we get 

K^0^ + 2)M. (26) 
^ 2 ' 16a ^ ^ 

ho can be arbitrary. Then by substituting (26) into (22), we obtain 

6hl 3ho(a + lf 2c 3K2(a + l)\^ 4. 

. (27) 

It is interesting to investigate the (in)stabihty of the scalar theory (16) 
with potential (27). The perturbations from the solution (19) are: 

H ^ho + hi{ts-t)" + Sh , (1)^ (poits-t)"^ +54) . (28) 
Then from the FRW equation (21) and 0-equation (24), it follows 
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Here 
Since 



Ho = ho + h{ts-tf , % = (l)o{ts-t)— . (31) 



^ l^oj = ^ (Is - tj 2 H (ts - t) 



a-3 
2 
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T/" ^ris ^ 3(q;- l)/io 1 3K:^(n + l)-(3a - l)c6g _i 
1/ l^oj - ^ (Is - H ^^^^ (Is - t) 



Hence, when t — t^, if « > —1, the third terms of V {^o) and V^" ($o) 
dominate and if a < —1, the second terms dominate. For the case of Barrow 
model, < a < 1. When t ~ t^, Eqs.(29) and (30) are: 

12, ^, (a + 1)60 , dbS 
~ -h^bh - ^ (ts -t) 2 ^ 



a - 1) 00 
4 



a-3 



{ts-t)—5(t) (33) 



(^50 3(a+l)0o,, 

i^-^)i^-Ait,_trH. (34) 



4 

By deleting 5/i from (33) and (34), we obtain 

0^gy- (°-^)(°-^) fa-r-^^. (35) 

Its solution is given by 

5(j)^(j)^{ts-t)'^ +(i>2{ts-t)'-^ . (36) 
Here 01 and 02 are constants. Hence, 

^^_ 6(« + l)(a-2) 
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In this order of the perturbation, 0i does not appear in 5h. Since, of course, 
2^ < when t ^ ts, the first term in 5(p (36) becomes large more rapidly 

than the unperturbative part $o = {^s ~ t) ^ (30), which tells that the 
solution (19) which describes the sudden singularity model in scalar-tensor 
theory is not stable. Already on the classical level, such instability may stop 
the appearence of future singularity. However, more secure mechanism which 
acts against the singularity occurence is quantum effects account. 

Near the singularity at t = t^, the curvature becomes large in general. As 
the quantum corrections usually contain the powers of the curvature ( higher 
derivative terms), the correction becomes important near the singularity. 
One may include the quantum effects by taking into account the conformal 
anomaly contribution as back-reaction near the singularity. The conformal 
anomaly T4 has the following form: 

TA^h(^F+ ^Di?^ + h'G + h"UR , (38) 

where F is the square of 4d Weyl tensor, G is Gauss-Bonnet invariant. In 
general, with N scalar, spinor, A^i vector fields, N2 (= or 1) gravitons 
and A^HD higher derivative conformal scalars, 6, h' and h" are given by 

AT + my 2 + l'2Ni + 611 Ar2 - 8A^HD 
~ 120(47r)2 

Ar+llA^i/2 + 62Ari + 1411Ar2-28ArHD 
^ - 360(4^0^ ■ ^ ^ 

As is seen 6 > and 6' < for the usual matter except the higher derivative 
conformal scalars. Notice that h" can be shifted by the finite renormalization 
of the local counterterm E?, so h" can be arbitrary (in ref. [10] it was chosen 
to be zero, for simplicity). In terms of the corresponding energy density pa 
and pressure pa, Ta is given by = — + ^Pa- Then by using the energy 
conservation law in FRW universe 

= ^ + 3HipA + PA) , (40) 



we may delete pa as 



r. = -4,.-l^|i, (41, 
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which gives the following expression for pA'- 



Pa 



dta^HTA 



= r / dta^H 



fdH\ 



fdH\ 



,dH 



-6 ( -6 + b" 



■'''[-dF) ^''''[-[-df) 

d^H ^,,d''H JdHV ,^,,odH 



■ (42) 



Moreover, 



PA 



1 dp A 

m dt 

Ta 1 f 4 



dta'HTA 



1 



-(f) 


+ 246'|- 






J y" dta^H 


-(f) 







d'^H ^ fdH\ 



,dH 



df 

2 



+ 4 — +12i/^— 



V (it / 



dt ) dt 
d^H ^ fdHV ,^,,odH 



■ (43) 



As in (19), one assumes 



H r^h'^, + h[(ts-tf , or a = aoe'^o*-^^*^-*) 



c'-3 



(44) 



We also consider the case t tg and keep only the first and the last terms in 
V{<j>) (27) 

^ - fc±i£0,^^^(f^i) . 



(45) 



These terms are dominant for Barrow model < a < 1. Then the consistent 
solution is given when a' > 2. From (18), it is seen the behavior of is not 
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so changed from the classical solution (19) 

ct>^ct>oits-tf . (46) 
Now the quantum corrected FRW equations are: 

Substituting (44) into (42), we obtain 

PA ~ -2Ab'h'o%-^^'o' J dte^^'o' = -eb'h'o" + pAo . (49) 

Here pAo is the integration constant which may be chosen to be zero since 
Pa ^ when b' ^ (classical limit). Substituting (46) and (49) with pAo = 
into (47), we obtain 

= -b'n^h'o^ - h'o^ + hi , (50) 

which can be solved as 



, _ i±v^r±iEZM > (51) 



if 



1 + Ab'K'^hl > , (52) 

which gives a non-trivial constraint since 6' < in general. In (51), the minus 
sign in ± corresponds to the classical case (19) in the limit of b' — > 0. In 
(48), when t tg, one finds 



^^2UJ "^^^^ ^ ■ ^ ^ 

First interesting case is that |6 + 6" does not vanish. Since 



PA -2[h + b'^^'^ 



a' -3 



~ 2^6 + b"^ a' {a' - 1) {a' - 2) h[ {t^ - t)' 



(54) 
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the classical term + 3i/^) in (48) can be neglected as this term be- 

haves as {ts — t) 



a'-l 



a 



. Due to p ~ 

= a + 2 , h[ 



-pA, one gets 

^ {a + l)0g 

8 (§6 + 6") {a + 2)a 



As 



Hr^hQ + h\ (ts - t) 



a+2 



(55) 



(56) 



the singularity at t — tg is moderated, that is, the exponent of the power 
oi ts — t becomes larger. When |6 + b" > 0, near t — ts, H decreases 
with time, that is, the universe is deccclcrating. On the other hand, when 
|6 + b" < 0, H increases with time, that is, the universe is accelerating. 
If we may replace [ts — t) with its absolute value \ts — t\, the deccelerating 
(accelerating) universe turns to accelerate (deccelerate) when t > ts- 

Another interesting situation corresponds to |6 + 6" = by properly 
choosing b". In this case 



PA ~ -4 (6 - 2b') ' ^ -4 (6 - 2b') K'a" [ts - tf^'^'-'^ 



(57) 



The choice consistent with (48) is 

, a + 1 



a 



4>l 



4 (6 - 2b') 



(58) 



Since a' — a = > 0, the singularity is moderated, compared with the 
classical case (19). 

One may also consider the case that the classical energy density p and the 
pressure p can be neglected since the quantum induced pA and pA become 
significally dominant. In this case, combining the first FRW equation and 
(42) 



dta^H 



fdH\ 



fdH\- 



126 — + 246' <^ - — + H'— + H 



-6( -6 + 6" 



d^H 
dt^ 



+ 7H 



\dt J 

d^H 



\dt J 



dt 



fdH\ 



dH 



df^ 



+ 4 — + 12H^— 



\dt J 



dt 



(59) 
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one has 



12 / ^ „ dH\ 



fdH\ 



fdH 



dt 



-6[^b + b" 



d^H 
dt^ 



+ 7H 



dt^ 



+ 4 



dH\ 



dH 



— 1 + 12H^ — 
dt ) dt 



Notice Eq.(60) is nothing but 



,(60) 



-.R = -Ta 



(61) 



Eq.(60) has a special solution, which gives a deSitter space with constant H. 
In fact, if H is assumed to be the constant, Eq.(62) reduces to 



24 



-^H^ = -2Ab'H^ , 



which has solutions 



1 



(62) 



(63) 



The second solution describes deSitter space. 

If the curvature becomes significally large, one may neglect the classical 
part, which is the l.h.s. of (60). Assuming 



H 



ho 
t ' 



(64) 



one arrives at the following algebraic equation: 

= 12/io |-3 06 + b"^ + (^-b -2b' + 9 06 + b")) ho 
+ (^_2b'-6(^^b+b"^^hl + 2b'hl^ . 



(65) 



Since the part inside { } of (65) is the third order polynomial, there is always 
a nontrivial solution for ho, at least, if ^b+b" does not vanish. If the obtained 
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ho is negative, the universe is shrinking but if we change the direction of the 
time hj T ^ ts — t, we may obtain a solution describing the expanding 
universe. Even if |6 + b" — 0, one gets non-trivial (non-vanishing) solution 
for hg: 

Since b" is arbitrary in principle, we may consider the case that the terms 
with b" become dominant. Then Eq.(60) reduces to 

12 dH\ ,,>fd^H ^,,cPH JdHV ^^^odH 

' (67) 

which can be written by using the scalar curvature R — % (2i?^ -|- ^) as 

2 „ [d'^R dR\ 

Eq.(68) has been found in i?^-gravity [16, 17] with the purpose to describe 
the inflation. Thus, like in ref.[10] (where b" was chosen to be zero) we come 
to the following picture. Near to future singularity, the quantum effects 
become dominant and they drive (most probably) the universe to deSitter 
space. Thus, final state of such universe is not the singularity. Rather, 
far in future the new inflation era (which is supported by quantum gravity 
effects[ll]) starts. 



3 Final state of dark energy universe 

In the same way the singularity avoidance in other models (of dark energy 
universe) may be considered. First of all, let us give the simple argument 
stressing that Big Rip should not occur. Working in adiabatic approximation, 
one supposes that H is almost constant and the time-derivatives of H can 
be neglected. Then since a oc e^*, using (42), we flnd 

Pa 246' i/^ . (69) 

The first quantum corrected FRW equation looks like 

6 

-H^^p + PA = p-2Ab'H^ , (70) 
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The above equation can be rewritten as 



= -246' + 



+ 



(71) 



Since h' and therefore -^-^ are negative, in order that H'^ has positive real 
solution, it follows the constraint for p 



P< 



(72) 



Thus, even if p includes the dark contribution from phantom, p has an upper 
bound. In other words, it does not grow infinitely with the time, which 
was the disaster for phantom cosmology. Equivalent upper bound may be 
suggested when one uses Hawking radiation from cosmological horizon (as it 
was communicated to us by P. Wang). Of course, near the Big Rip singularity, 
the time-derivatives of H should be taken into account in the consistent 
treatment of the sort presented in the previous section. 

Now we consider the Big Rip singularity [5] generated by the scalar field 
with the exponential potential: 



-1 J d^x^(^R-^d,<pdy-v . 



(73) 



When 7 < 0, the scalar is a phantom with w < —1. By solving the 0-equation 
of motion 



= -7 
and the first FRW equation 



dt 



(74) 



6 , 



when 



V{4>) = l^oe" 



20 



one gets a singular solution: 



= 00 In 



t.-t 



70^ 1- 



2^0 



(75) 
(76) 

(77) 
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which gives 



a 



ts-t 

U 



4 



(78) 



Here a is singular at t = if 7 < 0. General solution of above phantom 
system has been found in [11]. Even for the general solution, the behavior 
near t — ts is not qualitatively changed from that in (78). Hence, from the 
first look the Big Rip singularity seems to be inevitable. 

Near the Big Rip singularity, since a blows up, curvature becomes large 
as R (X \t — ts\ ■ Since the quantum correction contains powers and higher 
derivatives of the curvatures in general, the quantum correction becomes 
dominant. Hence, one can apply the same reasoning as in the previous sec- 
tion. With the account of the quantum correction (42), the corrected FRW 
equation has the following form: 



6 „2 7 
"^^ ^2 



K, 



dt 



+ V{(j)) + PA . 



(79) 



Let us assume 



H — ho + 5h , (j) = (poln 



t,-t 



+ S(I) . 



(80) 



and when t — > tg, Sh, 5(f) are much smaller than the first terms but ^ can 
be singular. Then 0-equation of motion (74) reduces to 



-7 



3hn 



{ts-ty ts-t J (Po{ts-t) 

+o{{ts-t)-^) , 



which gives 



Votl^-^-f, 5<i>^-l{ts-t) 



With |6 + 6" 7^ 0, one gets 



(81) 
(82) 



(83) 
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Substituting (82) and (83) into the quantum corrected FRW equation (79), 
we find 

''^:b + b"]^ + o{{t,-tr') , (84) 



= ^:^ + 6/^0 



and 



ts-t 



700 



(ts - t) In 



t.-t 



to 



Here t2 is a constant of tlie integration. Tlie scale factor a behaves as 



a = ao 



ts-t 


4(§b+6") 


t2 





(ts-t) _ho{ts-t)- 

6 



(85) 



(86) 



There appear logarithmic singularities in ^ but the singularity is mod- 
erated. Moreover, the universe might develop beyond t = tg- Thus, qTiantum 
effects prevent from the most singular universe. In case |6 + 6' = 0, the 
assumption (80) seems to be inconsistent. 

Another interesting dark energy model (which describes current accel- 
eration and even current dominance of dark energy) was proposed in [13], 
where the matter Lagrangian density (dark energy) is coupled with the scalar 
curvature: 

1 



(87) 



Here is matter-like Lagrangian density. The second term may be induced 
by quantum effects as some non-local effective action. By the variation over 
Qij^i,, the equation of motion follows: 







1 5S 



-9 Sg, 



1 



(88) 



Here the effective energy momentum tensor (EMT) T^,y is defined by 



6 



d'^x^/^Ld 



Let free massless scalar be a matter 



Ld = 



(89) 



(90) 
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The metric (2) is chosen. Assuming only depends on t (0 = <p{t)), the 
solution of scalar field equation is given by 

= qa-^R"^ . (91) 

2 

Here g is a constant of the integration. Hence R^Ld — ^ , which becomes 
dominant when R is small (large) compared with the Einstein term ^-R if 
a > -1 (a < -1). 

The accelerating solution of FRW equation exists [13] 

a+i / a + 1 
a — aot 3 -n = 



3t 

K^q^ (2a -I) {a- 1) 
3(« + l)"+^ (I (2a -1))' 



al ^ "'^ 'Z~ " '^2 ■ (92) 



Eq.(92) tells that the universe accelerates, that is, a > if a > 2. 

For the matter with the relation p = wp, where p is the pressure and p is 

2 

the energy density, from the usual FRW equation, one has a oc f3(™+iT. For 
a oc t^" it follows w = —1 + gl^, and the accelerating expansion (/iq > 1) of 
the universe occurs if it; < — |. For (92), one gets 

1 ~ _ 

w = . 93 

1 + a ^ ^ 

Then if a < —1, w < —1, i.e. an effective phantom. 

When a < —1, i.e., w < —1, the universe is shrinking in the solution 
(92). However, if one changes the direction of time ast ^ tg — t, the universe 
is expanding but has a Big Rip singularity at t — tg- Since near the singu- 
larity, the curvature becomes very large again, we may include the quantum 
correction (42) 

+ p + (94) 

If H behaves as in (92), after changing the direction of the time as t ^ tg — t 
in a, 

1 da a + 1 

adt 3{t-ts) ' ^ ^ 

the quantum correction of the energy density pA behaves as pA ~ (t — ts)~^, 
which becomes very large when t ^ tg- This shows that H cannot grow as in 
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(95). If H is not very large, p (94) becomes very small when a is large and 
can be neglected since p oc a^^. In such a situation, Eq.(95) reduces to (59). 
Hence, instead of future singularity, due to quantum effects the dark energy 
universe ends up in the deSitter phase (63). Thus, quantum effects resolve 
the sudden future singularity of dark energy universe. 



4 Thermodynamics and entropy bounds in 
the dark energy universe 

Thermodynamics of dark energy universe was discussed in ref.[8] where the 
appearence of negative entropies for models with equation of state param- 
eter less than —1 was demonstrated and entropy bounds were constructed. 
In the present section the entropy bounds near the Big Rip singularity are 
considered. The Hubble entropy Sh, Bekenstein entropy Sb, and Bekenstein- 
Hawking entropy Sbh are defined by 

^^-2G' ^^^-2G;^- 

Here G is a gravitational constant [k^ = 16nG) and V is the volume of the 
universe where for the universe with fiat spatial part, it is chosen 

V = Voa^ . (97) 

Near the Big Rip singularity, the dark energy dominates and the usual 
matter contribution may be neglected. Without quantum correction, a ~ 
(t, -t)3(^ and p ~ a~^(^+"') ~ {ts - 1)'^ in accord with (77) and (78). 
The entropies behave as 

2(l-3m) 1-m 4 

Sb ~ {ts - t) , Sh ~ (ts - t) »+i , Sbh ~ (ts - t) , (98) 

where the exponents are related by 

2(1 -3w) 1-w 4 

r < < r < , 99 

3{w+l) w + 1 3(w+l) ' ^ ^ 

when w < —1. Hence, all the entropies are singular at t = tg. Eq.(99) 
shows that the Bekenstein entropy Sb is most singular while the Bekenstein- 
Hawking entropy is less singular. 
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In order to estimate the entropy, we consider the thermodynamical model 
], where the free energy corresponding to matter with w is given by 



(t^V^ . (100) 



Here T is the temperature and V is the volume of the system. F is a function 
determined by the matter. The thermodynamical parameters are 



P 



' dV 



^-T^+ip' (t^v^ , 



S = -^ = -f(t^V^-^Tt^F'(t^V^ . (101) 



Here <S is an entropy. Since p behaves as p = po^"^*^^'''"''' ^ Po {^) i 
from the second equation (101) it follows 

PoVo^+'" oc -- [T^V^ F' [T^V^ . (102) 

Then T'^V should be a constant, which indicates that the entropy S in 
the third equation (101) is also a constant. Since Hubble, Bekenstein, and 
Bekenstein-Hawking entropies (98) diverge at the Big Rip singularity, the 
following entropy bound holds near the (classical) Big Rip singularity 

S <Sbh<Sh<Sb . (103) 

Here Eq.(99) was used. Then all the bounds are satisfied, which may be 
compared with the case of the brane-world dark energy model of ref. [18], 
where the Hubble entropy bound is not satisfied and the Bekenstein bound 
is often violated. 

The Hubble parameter in the expanding universe is given by if = z{w+i) ^ 
if w > -1 and H = w^j^ts - t if w < -1, ^ < when w > -I and ^ > 

S^w+l) * 'at at 

\iw < —1, which corresponds to the Big Rip singularity. Sb oc H'^Va oc H^a^ 
as found from the FRW equation. Then one may define the following quantity 
(which may indicate the future singularity occurence) 

S^^, (104) 
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where S oc since Sbh oc a^. Hence, if 5" decreases with time, there might 
occur the Big Rip singularity. 

Let us reconsider the above entropy bounds with the account of the quan- 
tum effects. Using the solution (80), (82), (85), (86), we find the entropies 
(98) behave as 

Then Sh and Sbh are finite and they may give meaningful entropy bound but 
Sb^s negative since 7 < and diverges. Hence, the Bekenstein bound S < Sb 
is violated. In (105), however, we have included only the classical part (75) 
in order to estimate Sb- With the account of the quantum correction pA (42), 
the singularity in Sb can be cancelled. Since p = pcb + Pa = ^H"^ = we 
find the following expression of the quantum corrected Bekenstein entropy 

S« = MM , (106) 

which is positive and finite. With S*^ (106) instead of Sb (105), all the 
entropies are finite. We should note 

S^-:^ Sh-:^ Sbh if /lo > — , 

ao 

S's^Sh^. Sbh if h <^ — . (107) 

ao 

The parameters ho and Oq are the values of the Hubble parameter and the 
size of universe at t = ts, which may be determined from the proper initial 
conditions. 

In [21], the (quantum corrected) entropy bounds have been discussed. In 
[21], the spatial part of the universe is sphere, where we have a relation [22] 

Si + {Sbh - Sb)' = Si , (108) 

even with the quantum correction. In case that the spatial part is fiat, 
Eq.(108) reduces to S^ — 2SbhSb- We should note that for Sh, Sbh (105) 
and S% (106), it holds 

S1^2SbhSI. (109) 
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Even for the classical case that all the entropies are singular, Sjj = 2SbhSb, 
which can be found from the FRW equation (75). One can rewrite the 
FRW equation (75) in the form Sjj — 2SbhSb by using the definition of the 
entropies (98). 

To conclude, it is shown that entropies near to classical singularity are 
singular as well. However, quantum corrected entropies are finite and give 
the well-defined entropy bounds. This is not surprised due to the fact that 
quantum effects help to escape of future singularity in the dark energy uni- 
verse. 

5 Dark energy universe with general equa- 
tion of state 

So far we concentrated mainly on the various aspects of the dark energy uni- 
verse with negative equation of state parameter which is less than —1. Nev- 
ertheless, the recent astrophysical data admit also the case of time-dependent 
equation of state parameter. Let us consider several examples of such dark 
energy cosmology and its late time behaviour. Note that several models of 
dark energy universe with time-dependent equation of state were discussed 
in [8, 23] (see also refs. therein). 

One starts from the general equation of state of the form 



instead of the equation p — wp with constant w. In (110), / can be an 
arbitrary function. Imagine that solving gravitational equations, we want to 
construct the cosmology with time-dependent t«, which describes the transi- 
tion from the deccelerating universe to the accelerating one. As an example, 
the following scale factor a{t) may be considered 



P = f{p) , 



(110) 



a = aoe t' 



(111) 



Here A and a are some constants. Hence, 

TT f ^ 

— = ah = a A + 
at \ 




(112) 
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In the case that A and a are positive, the universe is accelerating if 



t>to^ , (113) 

and deccelerating if 

t<to. (114) 

That is, the deccelerating universe turns into the accelerating one at t = to- 
Thus, if transition point occured about 5 billion years ago, the solution 
may approximately describe our universe. Note that to is positive when 
< a < 1. 

By using (6), one finds 

P=^(A + f)\ (115) 
p = _^('f!(^ + 5^+3A=) . (116) 



Eq.(115) can be solved as 



t=^^, (117) 



substituting (117) into (116), we obtain 



Hence, with (118) as the equation of state, we arrive at a solution (111), 
where the deccelerating universe turns into the accelerating one. 

General case (110) may be considered as well. Using the first FRW equa- 
tion (4) and the energy conservation law 



one gets 



0^^ + 3H{p + p) , (119) 



^ = F(p)^-kJ5(p + /(p)) . (120) 
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With a proper assumption about function /(p), we can find the t-dependence 
of p by solving Eq.(120). Using the obtained expression for p = p{t), one can 
also find the t-dependence oi p as p— f {p{t)). 

By combining (110) and (120), the pressure p can be expressed as 

Therefore if F{p) > (F(p) < 0), it follows w < -1 {w > -1). Wc now 
assume F(p) = at p = po; where po is a particular value of p. We further 
assume that when p ~ po, F{p) behaves as 

F(p)-Fo(p-po)" . (122) 

Here Fq is a constant and n is a positive odd integer. If n 7^ 1, by solving 
(120), one gets 

p ~ Po + {Fo (1 - n) (i - to)r^ , (123) 

and if n = 1, 

p ~ Po + Ce^°* . (124) 

Here to or C is a constant of the integration. Then p goes to po only at \t\ — > 
00, which may indicate that the region with w > —1 could be disconnected 
with the region w < — 1. Instead of a positive odd integer n, one may start 
from 

n=^^, (125) 
m 

with an integer m. Then as in (123) the time-dependent energy density looks 
like 

( F "1 ™ 

P~Po+|^(t-to)| . (126) 

Thus, the region w > —1 might be connected with the region w < —1. In 
this case, however, the equation of state has branches. 

We now consider the case with a linear equation of state p = wp where 
w depends on time as w — w{t). Replacing /(p) by w{t)p in (120) and using 
the first FRW equation (4), it follows 

2 / . 



H = -{ I {l + w(t))) . (127) 
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In order to investigate what happens when w changes the value from the one 
bigger than —1 to that less than —1, we now assume that near t = to w{t) 
behaves as 

w{t) ~ -l + woit-to) , (128) 

with constant Wq. Using 

j dt{l + w{t)) ~ ^wo {t - tof + wi , (129) 

one finds 

8 32 

P 7 9 r , H 5 . (130) 

3K^wo{t-to) +2wi} 3{wo{t-toy + 2wi} 

Here wi is a constant of the integration. Unlike the case in (123), there is no 
singularity aX t — to ii wi ^ 0. 

As one more example, the case with another w{t) may be considered 

= -1 - ■ (131) 

Here a, b, and to are positive constants. Then w{t) has the following proper- 
ties 

cit 

w{0) = wo = -l + -^> -1 , w{to) = -l, 
w{+oo) = Woo = -1 - a < -1 • (132) 
Hence, w{t) connects the region of to > —1 with that oi w < —1. Since 



W{t) = J dt {1 + w{t))^ a i^-t+{b + to) In^-^^ . 



(133) 



we find 



where ti is a constant of the integration. When t ^ to, one gets 

W{t) ^ a (^-to + {b + to) In ^) + + 0{{t- to)') , (135) 
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which is consistent with (129). If the universe is expanding, that is H > 0, 
at t = we find the following condition 

b>ti>0 . (136) 

W{t) behaves as W{t) ~ —at < when t is large. Thus, if the condition 
(136) is satisfied, W{t) vanishes at finite t { t — tg) where tg is a solution of 
the equation 

= W^(t,)=a|-t,+ (6 + io)ln^^| . (137) 

Hence, there appears singularity at t = > in p and H, which is 
nothing but the Big Rip singularity. However, even with (136), since 
— ^0 = (^ + ^o) In — h can be negative in general, the singularity may 
occur in the region w > —1. 

To conclude, we presented several examples of dark energy universe with 
time-dependent (negative) equation of state. The possibility to have natu- 
rally accelerated universe phase (sometimes, as a transition from deccelera- 
tion) is shown. It is interesting that when time-dependent equation of state 
parameter is negative (not only less but even bigger than —1) the finite time 
future singularity occurs as a final state of such universe. Nevertheless, in the 
same way as it was discussed in second and third sections one can show that 
quantum effects prevent the evolution to such final state ( eventually driving 
the universe to the inflationary era). 

6 Discussion 

In summary, we discussed several aspects of phantom thermodynamics and 
the final state of the phantom dark energy universe. Despite the absence of 
consistent phantom energy theory, some general results look quite promising. 
In particulary, it is shown that finite time Big Rip singularity remains to be 
deeply theoretical possibility in classical phantom theory. The account of 
quantum effects (when universe evolves to the singularity and when curvature 
invariants grow) is done. As a result, it is proved that quantum effects 
moderate the singularity or it even disappears completely. (Note also that 
stability analysis[ll], and gravitational perturbations account [19] indicates 
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that that perturbations act against the Big Rip occurence.) Hence, it is 
unhkely that the final state of phantom universe is Big Rip. Rather, the 
final state is the initial state on the same time, because infiationary era 
may start again in the future. The resolution of Big Rip singularity resolves 
also several related phenomena. For instance, entropy bounds which are 
divergent near Big Rip become well-defined after the quantum corrections 
are included. Similarly, escape of finite time singularity means that black 
holes mass evolution is less dramatic than it was predicted (masses do not 
vanish to zero). 

It is expected that soon precise observational cosmology data will give 
more stringent bounds for equation of state parameter. At the moment, it 
is still unclear if it will lie at quintessence, or at phantom region. Moreover, 
it is quite possible that the smiling universe hides a number of surprises for 
us. Nevertheless, phantom universe remains to be the theoretical possibility 
which is not explored yet and which deserves some attention. 
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A The entropy of phantom universe 

One more weird property of the dark energy universe with w less than —1 is 
the strange behavior of the entropy. In fact, it was pointed out in ref.[8] that 
entropy of such universe is negative. Another proposal came out in ref.[9] 
suggesting to consider phantom fluid as kind of cosmological quantum fliud 
(as nuclear spin model, for instance) where negative temperature is admitted. 
(Note that the idea of negative temperature in cosmological context was 
discussed first by Vanzo-Klemm [20]). In this case, the entropy may be 
positive. 

Let us describe the relation between the entropy and energy of such dark 
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energy universe when the temperature is negative. Starting from the model 
[8], instead of (100) we consider the following free energy: 

F^ = 7rF((7r)^y) . (i38) 

If the temperature T is positive, 7 = 1 and if it is negative, j — —1. Simple 
calculation gives the pressure p, the energy density p, and the entropy S: 

p = -(7T)^+- F' f(7T)- \/ 



p = --{-fT)'^^F'({^T)^v) , 
w V / 

F ((7r)» + ^T^F' ((7^)™ I . (139) 



w 
-7 



If the energy is extensive, the energy behaves as £^ = pV — > XE under the 
rescaling the entropy and the volume as <S — > XS and V — > XV. In accord 
with [8] we consider the following free energy: 

Fy, = -/o (iT)^^- V(l + fi (7r)-^ y-^) . (140) 

If there is no the second term, the first term gives the extensive energy. It 
is assumed the second term is small compared with the first term. Then one 
gets 



E = ?^ 

w 



k 

w 



fol{lT)^V{{l ^ 



w 



+ (l + -- — ]fi{^T)-^V-^] . (141) 
\ w nw J ) 

The sub-extensive part of the energy Eq., which is called the Casimir energy 
[22], is given by 

Ec = n{E^pV -TS)^-nV-'^(^ 

= -2/0/1 (7r)^+^-^ V'-^ . (142) 
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The extensive part of the energy Ee has the following form: 
Ee — E — -^Ec 



2 



= ^ (7r)^+- V (^1+ \^l-^ + wj fi {jT)-^ V-nj . (143) 

As in Section 4, T^V is a constant in the phantom dominated universe. 
Then if one neglects the second term in Ee and/or E as 

Ee^E^ ^(7r)^+- , (144) 

we obtain 

_ n 

(Kj+l)n-l 



S ^ A 



V'"^{2E- Ec)Ec . (145) 



Here A is a constant. 

The natural assumption is > 0. From the expression (141), /o < if 
w < 0. Let the starting condition is that the entropy S is positive. In case of 
the quintessence, where —l<w<\, since 1 + ^ < 0, from Eq.(141), we find 
7 > (positive temperature) so that the entropy is positive. On the other 
side, in case of phantom, where w < —1, that is, 1 + > 0, it follows 7 < 
if the entropy S is positive. Therefore the temperature should be negative. 
Conversely, if we assume the temperature is positive in the phantom theory, 
the entropy should be negative. 

Note also that in order to obtain Cardy-Verlinde formula (145) (for list of 
references, see review [24]), the Casimir energy Ec should be positive, which 
requires /i > 0. Hence, the entropy of phantom-filled universe is positive 
when the temperature is negative. In this case, standard CV entropy formula 
holds. 



B Black holes mass evolution in the dark en- 
ergy universe 

One more strange feature of the phantom universe is the black holes mass loss 
up to the full disappearence in the Big Rip singularity. The corresponding 
analysis [7] was performed in classical phantom-like universe (where domi- 
nant energy condition is broken) with final state in the Big Rip. In the present 
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Appendix, we reconsider this process taking into account the quantum effects 
which prevent the creation of Big Rip singularity as well as time-dependent 
(negative) equation of state. 

As is shown in an important paper [7] (see also [9]), the rate of the black 
hole mass change in the fluid with the energy density p and the pressure p is 
given by 

^ = 47rAM2 (p + p) . (146) 
at 

Here M is the mass of the black hole and A is a dimensionless positive 
constant. As a background, FRW universe with the metric (2) may be con- 
sidered. Combining the first FRW equation (3) and the energy conservation 
law (119), one obtains 

2 [2^(p^ 



Further combining (146) and (147), we get 



d f l\ %TiA 2 



d (^p2 

dt\Mj n V 3 dt 
The solution of above equation is: 



(148) 



M = ^^^^ . (149) 

I SttAMq /2p 

K y 3 

Hence, if p increases as in the case that the fluid is phantom, M decreases. 
At the Big Rip singularity where p diverges, M vanishes. This is an universal 
property for any black hole in such phantom universe. On the other hand, in 
the case of Barrow model where p is finite (5), the mass M is finite even at 
the singularity. By using the first FRW equation (3), we may further rewrite 
M (149) in the following form 

In [7, 9], the behavior of M for the phantom with constant w < —1 has 
been investigated in detail. As in (77), when w < —1 in the expanding 



30 



universe, the Hubble parameter behaves as 



H^--^^^^. (151) 
On the other hand, in case w > — 1 it looks hke 

2 

H = . (152) 

Using Eq.(150) one arrives at 

M = ^ — . (153) 

1 _ 4-irAMo 3{vj+i) 
3 ts-t 

when w < —1 and 

M = ^° , . (154) 

when w > — 1. In case of (153), M decreases and vanishes a,t t — tg. Near 
t — tg, M behaves as 

^^_9(.,+ l)fe-t) 

SttA ^ ^ 

This does not depend on Mq and is universal as pointed out in [7]. On the 
other hand, M (154) increases and reaches the maximal value M = Mq when 
t = oo. Even with the account of cosmological term, the qualitative behavior 
does not change. 

The above behavior is modified when quantum effects are taken into ac- 
count because as it was argued in second and third sections they may stop 
the evolution to final singularity. Indeed, let us consider the case that the 
quantum correction is included as in (79). If |6 + b" ^ 0, combining (80) 
and (85), it follows that H = ho when t = tg. Therefore M has a finite, 
non- vanishing value: 

Af 

^ - ■ <^«^' 

Since 

In 



dH d6h 700 
"dt ~ ~dt ~ ~2 {Ib + h") 



- 1 



+ 1 . (157) 
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and 7 < for the phantom, near the singularity ^ >0 (^<0) 

if |6 + b" < (§6 + b" > O). Therefore if §6 + b" < 0, since H increases, 
M decreases towards the singularity although M is finite and non- vanishing 
there. On the other hand, if |6 + b" < 0, M increases. As the Hawking 
radiation occurs due to the quantum correction, the above type of behavior 
may be more realistic in the phantom universe (quantum effects have been 
neglected, at least in the leading-order, in [7]). 

Let us reconsider what happens with the black hole mass M evolution 
equation (149) in the dark energy universe with time-dependent equation 
of state. For simplicity, the quantum corrections are neglected. When p — 
w{t)p, from (147) it follows 

(l+«,W)p=--^-^. (168) 

First we consider the case (128) and the behavior of p (130). Then if wq and 
wi arc positive, p takes a minimum value at t = to. The black hole mass M 
(149) increases when t < Iq and it reaches the maximum at t = to- When 
t > toj the mass decreases. 

As a more concrete example, w{t) (131) may be discussed. When ts > to 
(at t — ts, W{t) (137) vanishes), W{t) increases when t < to and decreases 
when t > to- Then the energy density p decreases when t < to, increases 
when t > to, and diverges at t = tg- Therefore the behavior of the black 
hole mass M (149) is similar to W{t), that is, M increases when t < to, 
decreases when t > to and vanishes at t = tg (like in classical phantom 
universe). Nevertheless, the account of quantum effects, as we showed above, 
qualitatively changes the black hole mass evolution. In other words, the same 
phenomenon which drives the dark energy universe out of final singularity 
(because of second quantum gravity era) is responsible for much less sharp 
loss of black holes masses. As Big Rip does not occur, initially massive black 
holes continue to be (may be less) massive! 
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